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INTRODUCTION 


Propellant sloshing is a potential source of disturbance that is critical to the stability of launch 
vehicles. The driving sloshing forces of the propellant masses come from the lateral 
disturbance, oscillatory thrust force, or angular rotation of the vehicle during the powered phase 
of the flight. Sloshing is expected to occur during vehicle taxi, takeoff, engine shut off, and 
flight maneuvers. Fluid motion produced during these operations can circulate sub-cooled 
liquid fuel near the liquid vapor interface resulting in increased condensation and corresponding 
pressure collapse. Conversely, liquid contact on hot tank walls caused by sloshing can result in 
rapid vaporization and subsequent rapid pressure rise. Understanding the dynamic behavior and 
the response of cryogenic propellants of Liquid Hydrogen (LH2) and Liquid Oxygen (LOX) 
tanks under spacecraft flight operations (such as engine restart, oscillatory thrust and oscillatory 
side forces) is extremely important during the vehicle design. In addition, the proper orientation 
of propellant is necessary for improved performance of the tank vent system, to ensure sufficient 
liquid propellant for engine firing and proper thermal conditions in feed lines prior to ignition, 
and to minimize vehicle disturbances and avoid excessive demands on the attitude controls. 

Sloshing is a wave phenomenon, and is driven by gravitational forces, surface tension, and 
surface kinematics. To take all the detailed dynamics into account in a vehicle dynamic 
modeling is very time consuming. Previously, a spring mass damper model has been used for 
Guidance, Navigation and Control Analysis (GN&C). Using first principle based Newtonian 
fluid dynamics; NASA Marshall Space Flight Center has developed a high fidelity 
Computational Fluid Dynamics (CFD) model for simulation and analysis of sloshing 
phenomenon. This model takes into account all the essential physics of fluid-gas interface, effect 
of gravitational force, surface tension and surface kinematics, to model the tank slosh using the 
acceleration forces provided by GN&C. 

The objective of this study is to validate CFD extracted mechanical parameters (such as sloshing 
mass, sloshing frequency, and sloshing mass center) against theoretical solution derived in 
1960’s for the simple cylindrical tanks and against experimental data taken during the 1960’s for 
more complex geometries 

2. ANALYSIS MODELING APPROACH 

2.1. Fluid Dynamics Analysis 

The basic CFD software used to study the tank sloshing phenomenon is the commercially 
available CFD-ACE+ code, which was originally developed by CFD Research Corp., and is 
currently owned and distributed by ESI Group USA. Some modifications and improvements, 
which will be described in the following section, have been made to tackle the unique vehicle 
sloshing problem. CFD-ACE+ has a VOF (Volume of Fluid) module which is designed for 
applications involving two immiscible fluids. A single set of momentum and continuity 
equations is solved, but different property sets are defined for each fluid. The volume fraction 
of one phase (in this case LH2 or LOX) is tracked throughout the solution to determine which 
fluid occupies each computational cell at any given time. In cells containing both fluids, a 
special routine is used to locate the shape and location of the interface. When there exists a 
surface tension force, its effect is applied in a conservative form. For time dependent 
simulations, such as unsteady tank sloshing, a special high order algorithm is used to update the 
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volume fraction in a cell from one time step to the next. A geometric reconstruction scheme for 
the interface representation is employed to track the interface. This unique algorithm is capable 
of capturing and maintaining a sharp interface between the phases for long time without 
introducing numerical diffusion. The time step size is determined by the local CFL number. 

It should be emphasized that accurate modeling of interfacial flows requires high fidelity 
algorithm for the kinematics and dynamics of interfaces, especially when surface tension force is 
important. In CFD-ACE+, the Laplace- Young high-order nonlinear boundary conditions at the 
liquid surface are applied, and a conservative algorithm for surface tension force previously 
developed by Yang and Przekwas [1] was implemented. This conservative surface tension 
treatment has shown to give high numerical stability, and can accurately predict transport speed 
of the liquid front, which is very difficult to do with other solution techniques. 

2.1. Implementation of Body Forces due to Non-Inertial Frame of a Vehicle 

The Navier-Stokes equations for a CFD code are typically written in the inertial reference frame 
and are derived from Newton’s law of motion in the form of: 


This law claims that the only way an object can change its velocity is if a non- zero net force is 
applied to the object. During the Crew Launch Vehicle (CLV) ascent, the Upper Stage (US) 
propellant tanks are in a non-inertial frame, in that the tanks are travelling along a curved path, 
and also speeding up or slowing down, and are subjected to continuous maneuvering including 
roll, pitch and yaw. In the non-inertial frame, Newton’s laws of motion of equation (1) will no 
longer hold, and there arise fictitious forces. This section will derive the extra fictitious forces 
in the non-inertial frame for the N-S equations. 


ma = f 


( 1 ) 


in: inertial frame, r: rotating frame 


inertial trame 

Figure 1. Inertial frame and non-inertial frame (rotating frame) 



In the non-inertial frame, any vector A has a relation of: 



( 2 ) 
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The subscript “in” refers to the inertial frame, and “r” refers to rotating frame (or non-inertial 
frame). As shown in Figure 1, 


Then 


let : r'= R + r, 


,dr\ .dR. .dr. _ - - 

^ =( ^ + ^ =v - +v - +axr 


( 3 ) 

( 4 ) 


Where, v cg denotes the translational velocity of the non-inertial frame origin (as observed in the 
fixed frame), v r is the velocity in the non-inertial frame relative to the origin of the frame. Now, 


, _ d ,dr\ _ ,dv 

let : a = — ( — ) ,a =( ) Jnen 

in 7 V 7 Jin* r V 

at at at 

dv cg x dv deb. „ _ .dr. 

a in ~ + (~ 7 ~)in + (~r)in X V + 03 X 

dt dt dt dt 


= a cg + (a r + a>xv r ) + aixr + a>x(v r + ®xr) 
= 5 cg +fl r +2fflxv r +fflxF + fflx(fflxf) 


( 5 ) 

( 6 ) 


It is interesting to note that: 

a cg : acceleration of vehicle center of gravity (CG) 

a r : acceleration of fluid particles in the non-inertial frame 

2 cbxv/. Coriolis acceleration 

cox r : angular acceleration 

cb x (cb x r ) : centrifugal acceleration 

With the above expressions, we can now write the Newton’s Law as: 

md in = f 

or: m {a cg + + 2<S xv r + axr + cbx(cbxr) = f 

ma r = f ~ma cg - mcbxr -2mcbxv r - mcbx(cbxr) 

FinaUy ma = f + f ( 7 ) 

r J J user 

Here f use r assembles the fictitious body forces due to vehicle acceleration and maneuvering, and 
they are the terms to be added into an existing CFD code. Specifically: 

(fuser ) X = ~ m ( a cg ) x ~ m(cb y r z - cb z r y ) - 2 m(a> y v„ - co z v ry ) - m[-(co y 2 + co z )r x 2 + co x o) y r y + co x co z r z ] 

(fuser )y = ~ m ( a cg )y ~ m (®J x ~ <*> x r z ) - 2 ™(o) z v rx - C0 x v rz ) - m[-(co z 2 + co 2 )r 2 + (0 y (0 z r z + co y co x r x ] ( 8 ) 

(fuser) z=~ m ( a c S )z ~ m (® x r y - cb y r x ) - 2 * 1 ( 0 , - co y v „ ) - m\-(co 2 + co 2 )r 2 + co z co x r x + co z co y r y \ 
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The above derived body forces due to non-inertial frame have been implemented into a user 
subroutine, which is then linked to the commercial CFD software, CFD-ACE+. 

3. VALIDATION STUDIES 


Sloshing is a wave phenomenon, and is driven by gravitational forces, surface tension, and 
surface kinematics. For LH2 and LOX propellant tanks, one can typically observe two types of 
surface sloshing waves: the first natural mode and the second natural mode. This is illustrated 

in Figure 3 for the LH2 tank and in Figure 4 for the LOX tank. 



First Mode, A=2D 



Second Mode, \=D 


Figure 3 The First Two Natural Modes Of LH2 Tank. D: Diameter Of The Tank, A: Wavelength, and g: 

gravitational vector. 





First Mode, A=2D 

Figure 4 The First Two Natural Modes Of LOX Tank. 


Second Mode, A=D 
D: Diameter Of The Tank, A: Wavelength. 


These two modes are standing modes and have a direct relationship with the diameter of the 
tank. The first mode is symmetric with respect to the central plane, while the second mode is 
axisymmetric. Typically, in flight, the first mode is the dominant mode. It is important to 
predict the natural frequency, the participating sloshing mass, the damping, the mass center of 
these modes. 


During the 1960’s when CFD technology was in its infancy, the sloshing dynamics of liquid tank 
was studied either analytically or experimentally. In the analytical approach, assumptions have 
been made to arrive the closed form solutions. While in the experimental approach, only sub- 
scale models are tested and measured. The full scale sloshing mechanical parameters have to 
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be extrapolated from the sub-scale model. To build confidence in using CFD technology to 
extract and analyze the sloshing dynamics, this paper will describe our effort to validate the CFD 
model against both analytical derivations and experimental measurements. 

3.1. Validation Against Analytical Solutions 

Southwest Research Institute has been involved in sloshing research since before the start of the 
US space program. Dr. Norm Abamson, then the leader of the sloshing group has organized 
and published a monograph [2] which contained very comprehensive knowledge on both 
experimental and analytical results of sloshing and related subjects. The analytical solutions 
were derived for simple geometries such as rectangular boxes or a cylindrical tank with a flat 
bottom face. These analytical solutions have been currently used at MSFC to get mechanical 
parameters (mass, frequency, damping, mass center) for the preliminary design and for sloshing 
model input into GN&C’s MAVERIC model. In reality, the analytical results are applicable only 
to the simple geometries such as cylindrical and spherical tanks. We will use the example of 
cylindrical tank to compare CFD predicted mechanical parameters with the analytical solutions. 

Simulation Model . The computational model along with the computational grid is shown in 
Figure 5. Here a cylindrical tank with diameter D=3” and height of H=ll” is modeled. The 
tank has no internal obstruction and all the walls are smooth. The top of the tank is open to 
ambient condition and is set up as fixed pressure boundary condition. All other boundaries are 
set as non-slip wall. There are a total of 80,000 cells, with 60 cells covering the whole tank 
diameter. This ensures that there will be at least 120 cells to resolve the first natural mode and 
60 cells to resolve the second mode. Butterfly grid is used to avoid the singularity at the tank 
center line, and to give better uniformity in the grid size. To determine the natural frequency of 
the tank slosh, an initial inclined surface is specified with a slope of 5.7 deg, as shown in Figure 
5. A typical time dependent mass center coordinate of the liquid mass inside the tank, 
computed by integrating the Volume of the Fluid value in each cell, is shown in Figure 6. 



Fixed P 


-slip 


Non-slip wall 



Cross-sectional butterfly grid 



Figure 5. CFD model of a cylinder tank sloshing for comparison with analytical solutions. 
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Figure 6. The center coordinates of the total liquid mass in the tank 

As seen from Figure 6, the vertical coordinate of the mass center is essentially constant, implying 
that the sloshing mass center does not move in the vertical direction, which is due to the 
symmetrical properties of surface wave. One may also notice an almost zero value of the mass 
center in the direction normal to the sloshing plane. This indicates that sloshing is essentially 
two dimensional. 

Extraction of Natural Frequency 

From Figure 6 one can see that the sloshing has a rather constant period, which makes it possible 
to extract the natural frequency of the tank sloshing from the above data. Several simulations 
were made at different heights of the liquid level, and they are summarized in Table 1. For the 
same configuration, one can use the analytical solution to predict the natural frequency as [3], 
which is also presented in Table 1: 





R 


R 


Table 1. Sloshing Frequency Comparison for a Cylinder Tank in Water, R=3” 


Liquid Level 

h/R=2.0 

h/R=1.0 

h/R=0.5 

h/R=0.25 

Analytical 

2.4523 Hz 

2.3880 Hz 

2.0867 Hz 

1.6058 Hz 

Present CFD 

2.4539 Hz 

2.3753 Hz 

2.0920 Hz 

1.6155 Hz 

Error 

0.06% 

0.53% 

0.56% 

0.97% 


Table 1 also shows the analytical solution, the CFD predictions and the difference presented in 
percentage error. One can clearly see that the CFD prediction gives excellent agreement with 
the analytical solution for all the liquid levels. This builds up the confidence of the CFD 
analysis. 
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Extraction of Slosh Mass and Slosh Mass Center from CFD Solution 

In addition to extracting the sloshing natural frequency, one can also extract sloshing mass and 
the sloshing mass center. The modeling procedure is to apply a periodic side load to force the 
liquid in the tank into a period motion at the forcing frequency. From the mass center 
amplitude, one can determine the equivalent sloshing mass of the system. The procedure and 
sample results are illustrated in Figure 7. Table 2 shows the comparison of the CFD prediction 
and analytical solution for the slosh mass ratio. Again, one can clearly see the results are in 
excellent agreement. 



Fs, de =agf*sin(27jfif) 


Extraction of Slosh Mass Procedure 

1 . Apply an oscillatory side force to the 
tank at frequency f 0 , and amplitude ag. 
other than natural frequency f n 

= ag*waffrf 0 () 

2. Once the mass center location reaches 
periodic state with maximum amplitude 
x max , determine the mass ratio from: 



Sample Mass Center History 


5 -(*£)’ 

m ag 


Figure 7. Procedure to extract sloshing mass. 


Table 2. Sloshing Mass Ratio (ml/m) Comparison for a Cylinder Tank in Water, R=3” 


Fiquid Fevel 

h/R=2.0 

h/R=1.0 

h/R=0.5 

h/R=0.25 

Analytical 

0.22735 

0.43236 

0.66028 

0.78252 

Present CFD 

0.22858 

0.43380 

0.66402 

0.77950 

Error 

0.54% 

0.33% 

0.56% 

0.39% 


3.2. Comparison with Experimental Data 
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The model is a 1/3.75 scale oblate spherical Centaur LOX tank, and experiments were carried 
out during the 1960’s by Sumner et al. [3], As shown in Figure 8, excellent agreement was found 
between the present CFD prediction and the experimental measurement. 


1/3.75 scale model of oblate spherical Centaur LOX tank 



Figure 8. Comparison of CFD extracted frequency vs. experiments for 1/3. 75 scale model of 

Centaur LOX Tank. 


3.3. Advantage in Using CFD Predicted Mechanical Model 

The previous section has clearly illustrated the soundness of the CFD approach in modeling 
detailed fluid dynamics of tank sloshing and with excellent accuracy in extracting mechanical 
properties for different configurations. This is especially applicable to the future design of 
propellant tanks, as there is currently no experimental data available for the same size and 
configuration. Table 3 lists a few points that CFD can help to improve the extraction of 
damped spring mass mechanical parameters. 
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Assumptions 
in Analysis 

Limitations of the 
Assumptions 

Advantages of CFD 
Modeling 

Tank Geometry 

• Flat bottom 
cylinders 

• Good for LH2 tank 

• Questionable for 
LOX tank 

• Realistically true tank 
geometry for Ares I tanks 

Liquid Height 

• Limiting shallow 
tank 

• Limiting deep tank 

• Good for LH2 tank 
during the first stage. 

• Not valid for LOX 
tank 

• No limitation 

• Applicable to all liquid 
heights during the first and 
second stages 

Sloshing 

Amplitude 

• Linear theory 

• Small amplitude 

• Large amplitude slosh 
could be encountered 

• Applicability to large 
amplitude slosh during 
flight is unknown. 

• Able to simulate the 
whole trajectory with small 
and large amplitudes. 

• Can shed light on the 
limitations of linear 
theory. 

Ullage Gas 

• No ullage gas or 
gas density is zero 

• No effect of ullage 
gas on the surface 
dynamics 

• Good for water and 
air. 

• Not good for 
cryogenic fluids with 
higher density ratio of 
gas to liquid. 

• Able to correct any error 
from the theory using the 
realistic properties, 
especially for the first 
sloshing frequency. 

Gravity (Body) 
Force 

• Constant 

• Neglected dynamic 
effect 

• g is a function of time 
during flight. 

• Dynamic effect are taken 
into account. 

• Can handle variable g. 


More detailed comparison with experimental data will be reported in the full paper. 
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Abstract 

Propellant slosh is a potential source of disturbance critical to the stability of space vehicle. The sloshing dynamics 
is typically represented by a mechanical model of spring mass damper. This mechanical model is then included in 
the equation of motion of the entire vehicle for Guidance, Navigation and Control analysis. The typical parameters 
required by the mechanical model include natural frequency of the sloshing, sloshing mass, sloshing mass center 
coordinates, and critical damping coefficient. During the 1960’s US space program, these parameters were either 
computed from analytical solution for simple geometry or by experimental testing for the sub-scaled configurations. 
The purpose of this work is to demonstrate the soundness of a CFD approach in modeling the detailed fluid 
dynamics of tank sloshing and the excellent accuracy in extracting mechanical properties for different tank 
configurations and at different fill levels. The validation studies included straight cylinder against analytical 
solution, and sub-scaled Centaur LOX and LH2 tanks with and without baffles against experimental results. This 
effort shows that CFD technology can provide accurate mechanical parameters for any tank configuration, and is 
especially valuable to the future design of propellant tanks, as there is no previous experimental data available for 
the same size and configuration. 


I. Introduction 

P ropellant sloshing is a potential source of disturbance critical to the stability of launch vehicles. The driving 
sloshing forces of the propellant mass come from lateral disturbance, oscillatory thrust force, or angular rotation 
of the vehicle during the powered phase of the flight. As such, sloshing is expected to occur during vehicle taxi, 
takeoff, engine shut off, and flight maneuvers. Fluid motion produced during these operations can circulate sub- 
cooled liquid fuel near the liquid vapor interface resulting in increased condensation and corresponding pressure 
collapse. Conversely, liquid contact on hot tank walls caused by sloshing can result in rapid vaporization and 
subsequent rapid pressure rise. Understanding the dynamic behavior and the response of cryogenic propellants of 
Liquid Hydrogen (LH2) and Liquid Oxygen (LOX) tanks under spacecraft flight operations (such as engine restart, 
oscillatory thrust and oscillatory side forces) is extremely important during the vehicle design. In addition, the 
proper orientation of propellant is necessary for improved performance of the tank vent system, to ensure sufficient 
liquid propellant for engine firing and proper thermal conditions in feed lines prior to ignition, and to minimize 
vehicle disturbances and avoid excessive demands on the attitude controls. 

Sloshing is a wave phenomenon, and is driven by gravitational forces, surface tension, and surface kinematics. 
To take all the detailed dynamics into account in a vehicle dynamic modeling is very time consuming. On the 
other hand, a complete dynamic-stability analysis of the vehicle must include the forces produced by the liquid 
motion in partially filled propellant tanks. The sloshing dynamics is typically represented by a mechanical model 
of spring mass damper. This mechanical model is then included in the equation of motion of the entire vehicle for 
Guidance, Navigation and Control Analysis (GN&C) analysis. The quantities that are input to the GN&C analysis 
to consider the propellant sloshing are: natural frequency of the sloshing, sloshing mass, sloshing mass center 
coordinates, and critical damping coefficient. Many analytical solutions have been made to determine these sloshing 
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mechanical parameters, but their accuracy and applicability to the real configuration are not well known for vehicle 
specific flight tanks. 

Using first principle based Newtonian fluid dynamics, NASA Marshall Space Flight Center has developed a high 
fidelity Computational Fluid Dynamics (CFD) model for simulation and analysis of sloshing phenomenon. This 
model takes into account all the essential physics of the liquid-gas interface, effect of gravitational forces, surface 
tension and surface kinematics, to model tank slosh using the acceleration forces provided by GN&C. The 
objective of this study is to validate CFD extracted mechanical parameters (such as sloshing mass, sloshing 
frequency, and sloshing mass center) against theoretical solutions derived in 1960’s for simple cylindrical tanks and 
against experimental data taken during the 1960’S for more complex geometries. 

II. Analysis Modeling Approach 


A. Fluid Dynamics Analysis 

The basic CFD software used to study the tank sloshing phenomenon is the commercial CFD-ACE+ code, which 
was originally developed by CFD Research Corp., and is currently owned and distributed by ESI Group USA. Some 
modification and improvement, which will be described in the following section, have been made to tackle the 
present unique vehicle sloshing problem. CFD-ACE+ has a VOF (Volume of Fluid) module which is designed for 
applications involving two immiscible fluids. A single set of momentum and continuity equations is solved, but 
different property sets are defined for each fluid. The volume fraction of one phase (in this case LH2 or LOX) is 
tracked throughout the solution to determine which fluid occupies each computational cell at any given time. In 
cells containing both fluids, a special routine is used to locate the shape and location of the interface. When there 
exists a surface tension force, its effect is applied in a conservative form. For time dependent simulations such as 
the present tank sloshing, a special high order algorithm is used to update the volume fraction in a cell from one time 
step to the next. The geometric reconstruction scheme for the interface representation is employed to track the 
interface. This unique algorithm is capable of capturing and maintaining a sharp interface between the phases for 
long time without introducing numerical diffusion. The time step size is determined by the local CFL number. 

It should be emphasized that accurate modeling of interfacial flows requires a high fidelity algorithm for the 
kinematics and dynamics of interfaces, especially when surface tension forces are important. In CFD-ACE+, the 
Laplace- Young high-order nonlinear boundary conditions at the liquid surface are applied, and a conservative 
algorithm for surface tension force previously developed by Yang and Przekwas [1] was implemented. This 
conservative surface tension treatment has shown to give high numerical stability, and can accurately predict 
transport speed of liquid front, which is very difficult to do with other solution techniques. 

B. Implementation of Body Forces due to Non-Inertial Frame of a Vehicle 

The Navier-Stokes equations for a CFD code are typically written in the inertial reference frame and are derived 
from Newton’s law of motion in the form of: 


ma = f (1) 

This law claims that the only way an object can change its velocity is if a non-zero net force is applied to the 
object. During the Crew Launch Vehicle (CLV) ascent, the Upper Stage (US) propellant tanks are in a non-inertial 
frame, in that the tanks are travelling along a curved path, and also speeding up or slowing down, and are subjected 
to continuous maneuvering of roll, pitch and yaw. In the non-inertial frame, Newton’s laws of motion of equation 
(1) will no longer hold, and there arise fictitious forces. This section will derive the extra fictitious forces in the 
non-inertial frame for the N-S equations. 
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in: inertial frame, r: rotating frame 



Figure 1. Inertial frame and non-inertial frame (rotating frame) 

In the non-inertial frame, any vector A has a relation of: 


dA dA - 


(2) 


With subscript “in” refers to the inertial frame, and “r” refers to rotating frame (or non-inertial frame). As 
shown in Figure 1, 


Then 


let : r'= R + r, 


,dr' 


,dR ,dr 




(3) 

(4) 


Where, v cg denotes the translational velocity of the non-inertial frame origin (as observed in the fixed frame), v , 
is the velocity in the non-inertial frame relative to the origin of the frame. Now, 


, _ d ,dr\ _ . dv 

let ’- a in = -rS~ 17 ) ^ a r =(—) r ,then 
at at at 

dv dv r . dd> . dr 

a. =( — — ). +( — -). +( ). xr+cux( — ). 

at at at at 

= a cg + (a r + co x v r ) + c5 x r + c5 x (v r + do x r ) 

= a cg + a r + 2(b x v r + co x r + d3 x (d3 x r ) 

It is interesting to note that: 

: acceleration of vehicle center of gravity (CG) 

a r : acceleration of fluid particles in the non-inertial frame 

2c5 x v : Coriolis acceleration 


(5) 


( 6 ) 


cox r : angular acceleration 

d5 x (c5 x r ) : centrifugal acceleration 

With the above expressions, we can now write the Newton’s Law as: 
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Finally 


(7) 


ma r = f - ma cg - md 5 x r - 2ma> xv r - mdd x (a> x r ) 

=f + fuser 

Here f user assembles the fictitious body forces due to vehicle acceleration and maneuvering, and they are the 
terms to be added into an existing CFD code. Specifically: 

0 fuser) x = ~ m ( a c g )x ~m(d) y r z - a z r y ) - 2m(a> y v a - co v ) - m[-(co y 2 +(o z 2 )r x +(o x co y r y +co x coj z ] 

2 2 ( 8 ) 

(fuser) y =-m{a cg ) y -m((0 z r x -ffl/ z )-2m(o z v„ -(o x v rz )-m[-((o z +co x )r y +co y co z r z +co y co x r x ] 
(fuser), =-m(a cg ) z - m(cb x r y - d> y r x ) - 2m(co x v ^ - (o y v J - m[-{(o x +co y 2 )r z +co z co x r x +co z co y r y ] 


The above derived body forces due to non-inertial frame have been implemented into a user subroutine, which is 
then linked to the commercial CFD software, CFD-ACE+. 

III. Validation Studies 

Sloshing is a wave phenomenon, and is driven by gravitational forces, surface tension, and surface kinematics. 
For LH2 and LOX propellant tanks, one can typically observe two types of surface sloshing waves: the first 

natural mode and the second natural mode. This is illustrated in Figure 2 for the LH2 tank and in Figure 3 for the 
LOX tank. 



First Mode, A=2D Second Mode, \=D 

Figure 2. The First Two Natural Modes Of LH2 Tank. D: Diameter Of The Tank, X: Wavelength, and g: 
gravitational vector. 
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First Mode, A=2D Second Mode, A=D 

Figure 3. The First Two Natural Modes Of LOX Tank. D: Diameter Of The Tank, A: Wavelength 

These two modes are standing modes and have a direct relationship with the diameter of the tank. The first 
mode is symmetric with respect to the central plane, while the second mode is axisymmetric. In the typical 
application, the first mode is the dominant mode. It is important to predict the natural frequency, the participating 
sloshing mass, the damping, and the mass center of these modes. 

During the 1960’s when CFD technology was in its infancy, the sloshing dynamics of liquid tanks was studied 
either analytically or experimentally. In the analytical approach, assumptions have been made such that a closed 
form solution can be obtained. While in the experimental approach, only sub-scale models are tested and 
measured. The full scale sloshing mechanical parameters have to be extrapolated from the sub-scale model. To 
build confidence in using CFD technology for the extraction and analysis of sloshing dynamics, this paper will 
describe our effort in validating the CFD model against both analytical derivations and experimental measurements. 


A. Validation Against Drop Tower Test 

First, the CFD-ACE+ code was validated for the sloshing problem of cryogenic propellant during spacecraft 
flight operation. The ground experiments were carried out at Marshall Space Flight Center using a drop tower 
facility to obtain a controlled low gravity environment for duration of up to 4.3 seconds [2]. The test was 
conducted in six-inch diameter cylinders (see Figure 4) using petroleum Ether as the test fluid. The test geometry 
in a planar presentation is shown in Figure 4(b). Due to symmetry, only half of the tank model is considered, as 
shown in Figure 4(c). 



(a). Testing Equipment (b) Cross-section. (c ) 3D model and grid 

Figure 4. Experimental Package, Test Model and Simulation Model 


The initial conditions, physical properties, and profile of reduction in acceleration with time are listed in Figure 
5. 
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Initial Conditions and Physical Properties 


Initial Conditions 


Physical Properties 



Density = 618 kg/m 3 
Surface Tension=0.0149 N/m 
Viscosity =0.0002488 kg/m-s 


Gravity Acceleration 

-9.8 m/s 2 0.0 <t< 0.08 s 

-1.75 m/s 2 0.08s <t< 0.09 s 

-1.06 m/s 2 0.09s < t < 0.10 s 
-0.89 m/s 2 0.10s < t < 0.20 s 
-0.51 m/s 2 0.20s <t< 0.30 s 
-0.21 m/s 2 0.30s <t< 0.356 s 
-0.10 m/s 2 0.356s <t< 4.0 s 


(0.0576m, 0.323m) 


Figure 5. Initial Conditions and Physical Properties of Modeled Drop Tower Test 


The computational results from the present CFD-ACE+ simulation, an experimental photograph, and a previous 
simulation using FLOW-3D [3] software are shown in Figure 6 at several selected time instances. One can see 
good agreement between the flow fields from CFD-ACE+ simulation and those of FLOW-3 D. The results also 
match well with the experimental data over the entire test duration. 


B. Validation against 2D Analytical Solution of Gravitational Wave 

In this section, an analytical solution is applied at the limiting condition of a 2D planar surface to validate the 
natural frequencies from the CFD simulations. Let us consider the two-dimensional gravitational wave shown in 
Figure. 7. 


Under the conditions of: 

1 . negligible viscous effect; 

2. 2D planar wave only; 

3. negligible gas density, or high density ratio between liquid and gas; and 

4. linear small perturbation wave, 

One can derive a dispersion relationship of [4]: 


co 2 = gk + 


ok 3 

P 


( 9 ) 


where oo is the angular frequency, g is the gravitational acceleration, k is the wave number, a is the surface tension, 
and p is the liquid density. With k=27i/^, X as the wavelength, we can find the frequency as: 


f = 


I g + 2no 
\ 2 ?d pX 3 


( 10 ) 


It is seen that the wave frequency depends on surface tension, gravity acceleration, and wavelength. A critical 
wavelength can be found when the surface tension contribution is the same as that of gravity: 
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CFD-ACE 




Experiment 










x 

Figure 6. Comparison of Experimental and Computational Flow Field Results 
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Figure 7. Surface Gravitational Wave on Liquid and Gas Interface 


This gives that: 


g 2 KG 
1k\ pX? 



(11) 

( 12 ) 


On the ground with LH2 properties, one can find that X c - 0.01m 
It indicates that surface tension is only important for very small ripples, 
surface tension contribution. This leads to: 



For this reason, we can neglect the 


(13) 


To verify that a CFD simulation can accurately predict the natural frequency, a simplified model was made as shown 
in Figure 8. 



First Standing Mode 

X=2*D~36 ft 


Second Standing Mode 

X=D~18ft 



D=18 ft 

Animation Animation 


Simplified 2D Model 

(can be used to compare 
with theoretical solution) 
Figure 8. Simplified two 
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Theory: f=0.3738 Theory: f=0.5287 

Prediction: f_water=0.3703 Prediction: f_water=0.5263 
f_LH 2 =0.3623 f_LH2=0.5 

dimensional model to verify the theory 
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First, the model is two dimensional planar, just as in the theoretical derivation. The simulation starts with an 
initial disturbed liquid surface profile either in the first mode or in the second mode (see Figure 8), and with the 
properties of water and air, due to their high density ratio. Once the waves settle down to their natural oscillation, 
the time required to complete a cycle is recorded, and the corresponding frequencies can be determined. The 
predicted first and second natural frequencies matched with the theory very well with errors less 1%. 

On the other hand, when using LH2 as liquid and gaseous hydrogen as gas, the natural frequency is slightly 
decreased in comparison to that of water and air as shown in Figure 8. This is very much expected, since the 
density ratio for LH2 is only about 10, and the gas density contributes the wave motion in the form of added virtual 
mass. This will make a lower frequency: from 0. 3703Hz to 0. 3623Hz for the first mode, and from 0. 5263Hz to 
0.5000Hz for the second mode. 

C. Sloshing Frequency Validation Against Analytical Solutions of 3D Cylinder Tank 

Southwest Research Institute has been involved in “sloshing” research since before the start of the US space 
program. Dr. Norm Abamson, then the leader of the sloshing group has organized and published a monograph [5] 
which contained very comprehensive knowledge on both experimental and analytical results of sloshing and related 
subjects. The analytical solutions were derived for simple geometries such as rectangular boxes or the cylindrical 
tank with flat bottom face. These analytical solutions have been widely used to get mechanical parameters (mass, 
frequency, damping, mass center) for the preliminary design and for sloshing model input into space vehicles. In 
reality, the analytical results are applicable only to the simple geometries such as cylindrical and spherical tanks. 
CFD can be become a powerful tool in predicting these parameters for realistic propellant tanks after systematic 
validation. In this section, we will use the cylindrical tank to compare CFD predicted mechanical parameters 
against the analytical solutions. 

Simulation Model 

The simulation model along with the computational grid is shown in Figure 9. Here a cylindrical tank with 
diameter D=3” and height of H=1 1” is modeled. The tank has no internal obstruction and all the walls are smooth. 
The top of the tank is open to ambient and is set up as fixed pressure boundary condition. All other boundaries are 
set as non-slip walls. There are a total of 80,000 cells, with 60 cells covering the whole tank diameter. This 
ensures that there will be at least 120 cells to resolve the first natural mode and 60 cells to resolve the second mode. 
Butterfly grid methodology is used to avoid the singularity at the tank center line, and to give better uniformity in 
the grid spacing. To determine the natural frequency of the tank slosh, an initial inclined surface is specified with a 
slope of 5.7 deg, as shown in Figure 9. A typical time dependent flow field and mass center coordinate of the 
liquid mass inside the tank, computed by integrating the Volume of the Fluid value in each cells, are shown in 
Figure 10. 




Cross-sectional butterfly grid 



Initial Condition 


Figure 9. CFD model of a cylinder tank sloshing for comparison with analytical solutions. 
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Total Mass Center Coordinate with Time 




Figure 10. The center coordinates of the total liquid mass in the tank and flow field. 



As seen from Figure 10, the vertical coordinate of the mass center is essentially constant, implying that the 
sloshing mass center does not move in the vertical direction, which is due to the symmetrical properties of surface 
wave. One may also notice almost zero value of the mass center in the direction normal to the sloshing plane. 
This indicates that sloshing is essentially two dimensional for this problem. 


Extraction of Natural Frequency 

From Figure 10 one can see that the sloshing has a rather constant period, which makes it possible to extract the 
natural frequency of the tank sloshing from the above data. Several simulations were made at different heights of 
the liquid level, and they are summarized in table 1. For the same configuration, the analytical solution of the 
natural frequency has been derived in [5], as 


/,= 


i 

In 


^ l-841g 


tanh( 


1.841/z 

R 


(14) 


Table 1. Sloshing Frequency Comparison for a Cylinder Tank in Wate r, R=3” 


Liquid Level 

h/R=2.0 

h/R=1.0 

h/R=0.5 

h/R=0.25 

Analytical 

2.4523 Hz 

2.3880 Hz 

2.0867 Hz 

1.6058 Hz 

Present CFD 

2.4539 Hz 

2.3753 Hz 

2.0920 Hz 

1.6155 Hz 

Error 

0.06% 

0.53% 

0.56% 

0.97% 


The resulting frequencies at different liquid fill levels calculated from the above relationship are also listed in 
Table 1. It is seen that the present CFD predictions show excellent agreement with the analytical solution for all the 
liquid levels. All the discrepancies are less than 1%. This builds the confidence of the CFD tool in predicting 
sloshing frequency. 

D. Sloshing Frequency Comparison with Experimental Data for Scaled Centaur Oxygen Tank 

An experimental investigation of liquid sloshing in a 1/3.75 scale-model Centaur stage liquid oxygen tank was 
conducted at NASA Lewis Research Center to determine (1) fundamental frequencies of oscillation, (2) the 
maximum horizontal slosh forces and damping ratios occurring at the fundamental frequencies, and (3) the slosh 
damping effectiveness of the proposed slosh baffle. Experimental data was obtained over a range of liquid depth 
ratios for the following tank configurations: (1) clean tank with no internal hardware, (2) unbaffled tank that 

included simulated internal flight hardware such as the trust distribution cylinder, and spring ring, and supporting 
structure in addition to the same internal hardware used in the unbaffled tank. The contained liquid was water in 
most tests. The results of this investigation were presented in terms of dimensionless parameters. 
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The experimental test facility is shown in Figure 11, and the CFD simulation model is also shown in Figure 11. 
The three dimensional CFD model has 275,000 butterfly cells to remove the singularity near the tank center. The 
1/3.75 scale model has a major diameter of axis of 32 inch, and minor diameter of 22.86 inch. A cylinder section 
with 0.88 inch height is located at the equator to simulate the cylindrical portion of the full scale tank. 
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J, Stitiodtry tutfcd 
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Test Facility CFD Simulation Model 

Figure 11. Experimental test facility and CFD simulation model for the Centaur liquid oxygen tank. 


Natural Frequency Comparison for Smooth Tank 

During the test, the tank was forced to oscillate sinusoidally at a preselected excitation frequency and amplitude. 
At each liquid depth ratio investigated, the excitation frequency chosen was equal to the fundamental frequency of 
the oscillation of the contained liquid to provide the maximum slosh forces at a given excitation amplitude. When 
the first mode waveform has built up to its maximum height on the tank wall, the oscillatory motion of the tank was 
quick-stooped, and the residual horizontal slosh forces were recorded on an oscillograph trace. The fundamental 
frequencies of the liquid oscillation were determined by using the first several slosh force peaks occurring 
immediately after the quick stop, and was put into the form of the fundamental frequency parameter: 

(i5 » 

where r is the radius of liquid surface of partially filled oblate spheroidal tank, and g is the gravitational 
acceleration (9.8 m/s 2 ). 

Both the experimentally determined frequency parameter and those from CFD modeling are given in Figure 12. 
As one can see, excellent agreements were found between the present CFD prediction and the experimental 
measurement. 
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Figure 12. Comparison of CFD extracted frequency vs. experiments for 1/3.75 scale model of Centaur LOX 
Tank. 

E. Sloshing Frequency Comparison with Experimental Data for Scaled Centaur Hydrogen Tank 

Centaur hydrogen tank is different from the oxygen tank in that it has a cylindrical top and an oblate spheroidal 
common bulkhead with the oxygen tank. The geometry of the tank along with the test facility is shown in Figure 13. 
Sumner et al [7] measured the fundamental frequencies of the liquid oscillations, sloshing force and damping ratio 
for both the unbaffled and the baffled tank configurations for a range of liquid levels above the top of the concave 
oblate spheroidal tank bottom. This configuration is unique in that no analytical solution can be derived for the exact 
same geometry. Tests had been done for the case with a baffle located at liquid depth ratio of 0.59, and this is a good 
validation case for CFD model in capturing baffle impact. 



A, 10 000-Pound preloaded horizontal load cell 

B, Preloading bracket 

C, 500-Pound verical load cells 

D, Support members from overhead crossbeams to frame 

E, 32.7-lnch-diameter cylindrical tank 

F, Testbed 

G, Frame 

H, Hydraulic cylinder 
X Stationary support 



Figure 13. Schematic of scaled-model Centaur liquid-hydrogen tank and test facility 


The pictorial view of the hydrogen tank and its model using CFD tool are shown in Figure 14. A baffle of 
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2.19” wide is located at liquid depth ratio of 0.59. The tank is discretized with a total of 410,000 cells. 



Experimental tank 

Figure 14. Test hydrogen tank and its CFD model 


CFD model and grid 


Natural Frequency for Smooth Tank 

The test procedure to obtain the fundamental frequency of slosh was similar to that for the oxygen tank. First, the 
tank was oscillated sinusoidally at a preselected excitation frequency (from 0.3 Hz to 4.8Hz) and amplitude. At each 
liquid depth ratio (see Figure 14), the excitation frequency chosen was equal to the actual fundamental frequency of 
the slosh to provide the maximum slosh force. When the waveform had built to its maximum height on the tank 
wall, the oscillation of the tank was quick-stopped, and the residual horizontal slosh forces were recorded on the 
oscillograph trace. The first natural frequency of the liquid oscillation was then determined by using the first several 
slosh-force peaks occurring immediately after the quick stop. The test data, CFD extracted frequency, and the 
analytical solution based on a straight cylinder are shown in Figure 15. The analytical solution based on the straight 
cylinder solution over-predicts the frequency parameter for the liquid depth ratio from 0.45 to 0.8. When the liquid 
depth ratio is above 0.8, the curved bottom effect can be neglected and a straight cylinder solution can be applied. 
On the other hand, the present CFD simulation can model different geometrical configurations of the tank, and 
shows very favorable comparison with the test data as seen from Figure 15. This validation case illustrates the 
capability of CFD in extracting mechanical parameters for complicated tank geometries when there is no analytical 
solution available. 

Natural Frequency for Hydrogen Tank with Baffle 

Next we consider the case when there is a baffle inside the tank. The baffle has a width of 2.9” and is located at 
liquid depth ratio of 0.59. Figure 16 illustrates the sloshing interface dynamics during CFD simulations at several 
liquid depth ratios. The comparison of natural frequencies from the present CFD and test data is given in Figure 17. 
The experimental test showed that the presence of the baffle sharply increased the natural frequency parameter when 
the liquid-depth ratio corresponded to the baffle location, and the baffle reduced the frequency parameter when the 
liquid level was slightly above the baffle location. This increase in frequency can be understood as the reduction of 
cross-sectional radius of the tank due to the existence of the baffle when the liquid level is at the baffle location. In 
general CFD simulation shows fair agreement with the test, including the sudden rise and fall observed in the test 
near the baffle location. Even though the predicted amplitude of rise and fall is not as dramatic as the test data, the 
fundamental physics of the baffle impact on the natural frequency is well captured in the present modeling. 
Systematic study is under way to investigate the discrepancy. 
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Figure 15. Comparison of sloshing natural frequencies for the scale-model of Centaur hydrogen tank 



Figure 16. CFD free interface dynamics at different liquid depth ratios inside scale model of Centaur 
hydrogen tank. 
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Sloshing Frequency Comparison 



Figure 17. Comparison of sloshing natural frequencies for the scale-model of Centaur hydrogen tank 
with a baffle 


F. Sloshing Mass Extraction Validation 

The above sections describe the effort in extracting fundamental frequency of the slosh from a CFD 
simulation. All results show good agreement in comparison with experimental data. For the mechanical model, 
one also needs the slosh mass and its center. In the experiment of reference [6], the tank was forced to oscillate at 
different amplitude to derive the wave amplitude of the liquid surface as high as possible. Once the wave height had 
reached a steady-state value, the oscillatory motion of the tank was quick-stopped, and the horizontal slosh force 
was determined from the first force peak occurring immediately thereafter to calculate the slosh mass. The similar 
procedure that can be used in CFD modeling to determine the slosh mass is by forcing the tank to a periodic motion 
with amplitude of body acceleration ag: 


ill 


s 



+ k s x s = -m s agsm(2noj f t) 


(16) 


where C0f is the forcing frequency. Under the steady state condition, the slosh mass will oscillate at a periodic 
state with a frequency the same as the forcing frequency C 0 f . Let: 


sin( 271(0, 0 


Noting that: 


m 


x t = 


-x 0 


m , 



(17) 


( 18 ) 


we can find that: 
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ag 


( 19 ) 


Using equation (18), we can find that: 



m, 


$o_ 

ag 


(o ) n 2 -a>/) 


( 20 ) 


Here, x t0 is the maximum mass center displacement during the cycle. Using this method, there is no need to 
compute the force acting on the tank walls. As shown in Figure 18, first a sinusoidal body force is applied to the 
cylinder tank of Figure 9. The liquid fill level is at h/R=2.0, and the natural frequency predicted from the above 
simulation and analytical solution is 2.45Hz. The forcing frequency is set at 1Hz, and the total mass center history 
in Figure 18 shows an initial transient and then a periodic motion matching with a forcing frequency. The maximum 
value of the mass center displacement is the value used for equation (20). 


Figure 18. Procedure to extract slosh mass and the sample mass center history 


Validation of Slosh Mass against Analytical Solution 

For the cylinder tank, the analytical solution of the sloshing mass to total tank mass ratio has been derived in [5], 


m, 


0.45477? t 1.84 h 
-tanh- 


h 


R 


( 21 ) 


Table 2 shows the comparison of CFD prediction and analytical solution on the slosh mass ratio, 
can clearly see the excellent agreement with error less than 0.6%. 


Again, one 
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Table 2. Sloshi ng Mass Ratio 

Liquid Level 


ml/m) Comparison for a Cylinder Tank i n Water, R=3” 


h/R=2.0 


h/R=1.0 


h/R=0.5 


h/R=0.25 


Analytical 
Present CFD 


Error 


0.22735 

0.22858 

0.54% 


0.43236 

0.43380 

0.33% 


0.66028 

0.66402 

0.56% 


0.78252 

0.77950 

0.39% 


IV. Conclusion 

The present study represented a continuing effort in validating CFD technology in modeling space vehicle 
sloshing dynamics [8,9]. This effort has clearly illustrated the soundness of the CFD approach in modeling the 
detailed fluid dynamics of tank sloshing and the excellent accuracy in extracting mechanical properties for different 
tank configurations and at different fill levels. The present study shows that CFD technology is especially 
applicable to the future design of propellant tanks, as there is no previous experimental data available for the same 
size and configuration. From the geometrical point of view, the analytical solution is only available for the straight 
cylinder. It is good fit for the LH2 tank when the liquid level is high. CFD modeling has no limitation, and can 
take the real geometry into account. The analytical solutions are generally obtained for the linear analysis, where 
the amplitude of the slosh wave is small. CFD solution can go beyond the linear stability regime, and predict what 
will happen after the surface breakup. CFD simulation has also no limitation on the gravitational acceleration, and 
can be applied to the realistic flight modeling. 

Extracting the realistic damping ratio is a much challenge task as the fluid damping depends greatly on the grid 
resolution, numerical scheme and the ability in identifying the boundary thickness a prior. This is especially true 
when there exists a baffle. The results of the validation study will be reported in the future. 
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